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ABSTRACT 



Accretion flows having low angular momentum and low viscosity can have 
standing shock waves. These shocks arise due to the presence of multiple sonic 
points in the flow. We study the region of the parameter space in which mul- 
tiple sonic points occur in viscous flows in the absence of cooling. We also 
separate the parameter space in regions allowing steady shocks and oscillat- 
ing shocks. We quantify the nature of two critical viscosities which separate 
the flow topologies. A post-shock region being hotter, it emits harder X-rays 
and oscillating shocks cause oscillating X-ray intensities giving rise to quasi- 
periodic oscillations. We show that with the increase in viscosity parameter, 
the shock always moves closer to the black hole. This implies an enhancement 
of the quasi-periodic oscillation frequency as viscosity is increased. 



MNRAS (in press) 



1 INTRODUCTION 

In the standard theory of thin accretion flows around black holes (Shakura & Sunyaev, 
1973, hereafter referred to as SS73) viscosity plays a major role. Viscosity transports angular 
momentum outwards and allows matter to sink into the potential well formed by the central 
compact object. In this model, the flow angular momentum is assumed to be Keplerian 
and this is the standard notion about how matter is accreted. However, Chakrabarti & 
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Molteni (1995, hereafter referred to as Paper I), and Lanzafame, Molteni & Chakrabarti 
(1998, hereafter referred to as Paper II), through extensive numerical simulations showed 
that the angular momentum distribution depends strictly on the viscosity parameter and 
the way the viscous stress is defined. They showed that close to a black hole, the disk 
does not have a Keplerian distribution. This is because the flow must be supersonic on the 
horizon (Chakrabarti, 1990a) whereas a Keplerian disk is always subsonic (SS73). In Papers 
I and II, it was shown that for a large region of the parameter space, shocks may form in 
accretion flows and when viscosity is increased beyond a critical value (Chakrabarti, 1990ab; 
Chakrabarti 1996, hereafter C96a), the shocks disappear. 

Paper I also improved the concept of viscosity parameter a (SS73): it argued that in 
a generalized flow with significant radial velocity the viscous stress w.^^ should not be 
equated to —aP as in SS73, where P is the total pressure, but to —au{P + p^^), (actually, 
its vertically integrated value using a thin disk approximation) where, p is the density and 
a subscript 11 is given to a to distinguish it from the Shakura-Sunyaev viscosity parameter. 
The latter prescription naturally goes over to the original prescription when radial velocity 
is unimportant as in the case of a standard Keplerian disk model (SS73), however, when the 
radial velocity is important as in the transonic flow solutions (Chakrabarti 1990a), the latter 
definition preserves the angular momentum even across axisymmetric discontinuities, such 
as accretion shocks. The reason is, according to the Rankine-Hugoniot conditions (Landau 
& Lifshitz, 1959), in a steady flow, the sum of thermal pressure and ram pressure, i.e., 
P + pd'^ is continuous across discontinuities. This makes the viscous stress Wr^, continuous 
across axisymmetric discontinuities as well. 

In an earher study, Chakrabarti (1989a, hereafter C89a) considered the transonic prop- 
erties of isothermal accretion flows and showed that for a large region of the parameter space 
spanned by the specific angular momentum and the temperature of the fiow, an accretion 
disk can have standing shock waves. The specific angular momentum of the disk was smaller 
than that of a Keplerian disk everywhere. This flows come about especially when the mat- 
ter is accreted from the winds of a binary companion. Subsequently, Chakrabarti (1990b, 
hereafter C90b) showed that inclusion of viscosity reduces the region of the parameter space 
in that, at a sufficiently high viscosity, the Rankine-Hugoniot conditions which must be 
satisfied at a steady shock are not satisfied anywhere in the fiow. Existence of standing 
shocks in sub-Keplerian inviscid accretion disks have been tested independently by several 
groups since then (Nobuta and Hanawa 1994; Yang and Kafatos, 1995; Lu and Yuan, 1997). 
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Properties of Accretion Shock Waves in Viscous Flows Around Black Holes 3 
Numerical simulations have also been carried out with several independent codes such as 
Smootherd Particle Hydrodynamics (SPH) and Total Variation Diminishing (TVD) and dis- 
tinct standing shocks were found exactly at the predicted locations (Chakrabarti & Molteni, 
1993; Molteni, Ryu & Chakrabarti, 1996). 

In more recent years, it has become evident that the standing shocks may be very impor- 
tant in explaining the spectral properties of black hole candidates (Chakrabarti & Titarchuk, 
1995, hereafter CT95) as the post-shock region behaves as the boundary layer where accret- 
ing matter dissipates its thermal energy and generates hard X-ray by inverse Comptoniza- 
tion. C96a considered unification of solutions of winds and accretion around compact objects. 
However, the cooling was treated in terms of a parameter and no parameter space was stud- 
ied. The post-shock region is also found to be responsible to produce relativistic outflows 
(Chakrabarti, 1999; Chattopadhyay and Chakrabarti, 2002). Furthermore, numerical sim- 
ulations indicated that the shocks may be oscillating at nearby regions of the parameter 
space in presence of cooling effects (Molteni, Sponholz & Chakrabarti, 1996) and the shock 
oscillations correctly explain intricate properties of quasi-periodic oscillations (Chakrabarti 
& Manickam, 2000). Recent observations do support the presence of sub-Keplerian flows in 
accretion disks (Smith et al. 2001; Smith, Heindl and Swank, 2002). 

In view of the importance of the sub-Keplerian flows we plan to re-investigate the work 
done on isothermal flow by C89a and C90b by extending them to study polytropic flows to 
check the properties of shock waves in viscous flows. What is more, unlike C89a and C90b, 
we investigate the behaviour of the solutions in the entire parameter space spanned by the 
specific energy, angular momentum and the viscosity. In C96a some work was done, the 
parameter space was not explored. We find very important results: even when the viscosity 
parameter is very high, the fiow continues to have three sonic points: a prime condition to 
have a standing or oscillating shock waves. However, the parameter space for standing shock 
waves is gradually reduced with the increase of viscosity. On the other hand, we discover 
that the shock location itself is reduced with the increase in viscosity parameter. 

We wish to emphasize that the problem at hand is by no means a trivial extension of 
previous works. In an accretion fiow, where the fiow is subsonic at a large distance and is 
necessarily supersonic on the horizon, the fiow has to first become supersonic at a sonic 
point and then after the shock transition where it becomes sub-sonic, the fiow must again 
pass through the inner sonic point before entering into the black hole. In studying flows with 
constant energy (Chakrabarti, 1989b, hereafter C89b) or an isothermal flows (C89a), both 
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the sonic points were known when the so-called 'eigen-values', namely, the specific energy 
(for polytropic flow) or temperature (isothermal flow) and the specific angular momentum, 
are supplied. In the present situation, neither of these two quantities is constant in the flow 
since the viscosity will heat up the gas, increase thermal energy and at the same time reduce 
the specific angular momentum as the flow proceeds towards the black hole. Thus, the inner 
sonic point, through which the flow will pass after the shock, is not known before the entire 
problem is actually solved. We have devised a novel way to solve the entire problem by 
iterating the location of the inner sonic point till the shock condition is satisfied. We have 
identified the topologies which are essential for shock formation. We have also identified the 
parameter space which will have solutions with three sonic points but need not have standing 
shocks. These solutions generally produce oscillating shocks as shown by Ryu, Chakrabarti 
and Molteni, (1997). In C96a, some effects of viscous heating were studied and cooling 
effect was chosen to be proportional to the heating effect for simplicity. No parameter space 
study was made. In the present paper, we ignore cooling completely. Exact effects of various 
cooling processes and their infiuence on the parameter space will be discussed elsewhere 
(Das & Chakrabarti, 2003). 

The plan of the present paper is the following: in the next Section, we present the model 
equations. In §3, we present the sonic point analysis. In §4, we study the global solution 
topology. In §5, we classify the parameter space in terms of whether a global solution has 
triple sonic points or not. In §6, we classify the region with triple sonic points further to 
indicate which region may allow standing shocks and which region may allow oscillating 
shocks in presence of viscosity. We showed in particular that matter with a very low angular 
momentum may allow shocks even when the viscosity parameter is very high. In C96a it 
was shown that topologies are changed with viscosity and there exists two critical viscosity 
parameters at which such changes take place. In §7, we quantify these critical viscosity 
parameters. Finally, in §8, we discuss about the relevance of shock shocks in the context of 
quasi-periodic oscillations and make concluding remarks. 

2 MODEL EQUATIONS 

We consider a steady, thin, viscous, axisymmetric accretion fiow on to a Schwarzschild 
black hole. The space-time geometry around a Schwarzschild black hole is described by 
the pseudo- Newtonian potential introduced by Paczynski & Wiita (1980). Here, one uses 
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the pseudo-Newtonian potential given by, g{r) = ~ ,,._2gmbh/c^ ■ consider the units of 
velocity, distance and time to be c, = 2GMbh/c^ and 2GMbh/c^ respectively where c 
is the velocity of light, G and Mbh are the gravitational constant and mass of the black 
hole respectively. In this unit, defining x — r/rg, we get the potential as g{x) — — 2(x-i) - 
assume the disk to be in hydrostatic equilibrium in the vertical direction. 

In the steady state, the dimensionless hydrodynamic equations that govern the infalling 
matter are the followings (C96a). 

(a) Radial Momentum equation : 

^(M ^ IdP _ Xjxf ^ 1 _ Q 
dx p dx x^ 2{x — 1)^ 

(b) Baryon number conservation equation : 



M = mx (16) 

apart from a geometric constant. 

(c) Angular momentum conservation equation : 

^d\ix) Id/ n^^^ \ 

ax L,x ax ^ ' 

and finally, 

(d) The entropy generation equation : 

Ei?r^ = Q+ - Q-. [Id) 

The local variables i?, p, P and \{x) in the above equations are the radial velocity, density, 
isotropic pressure and specific angular momentum of the fiow respectively. Here S and W^^ 
are the vertically integrated density (Matsumoto et al. 1984) and the viscous stress, s is the 
entropy density of the fiow, T is the local temperature. and Q~ are the heat gained and 
lost by the fiow (integrated in the vertical direction) respectively. 

In our model of the disk which is assumed to be in hydrostatic equilibrium in the vertical 
direction, local disk height is obtained by equating the pressure gradient force in the vertical 
direction with the component of the gravitational force in that direction. The half thickness 
of the disk is obtained as: 

h = ax^'^{x - 1). (2) 
Here, a is the adiabatic sound speed defined as a = \p)Pjp. As discussed in the introduction, 
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we shall use the viscosity prescription of Paper I valid rigorously for flows with significant 
radial motion. Thus the viscous stress is: 



where, 11 = W + 10}'^. As mentioned before, this will ensure that the viscous stress is contin- 
uous across the axisymmetric shock wave that we are studying here. It is to be noted that in 
SS73 prescription, Wx,/, is not continuous across the shock. Thus, the stress would transport 
angular momentum at different rates on two sides of the shock which would always move 
the shock one way or the other. This is unphysical, since in the absence of viscosity, a stand- 
ing, axisymmetric shock is perfectly stable. It is impossible that an infinitesimal viscosity 
should destabilize the shock. However, this would have been the case if SS73 prescription 
were rigorously correct. 



At the outer edge of the accretion disk, matter has almost zero radial velocity even though it 
enters into the black hole with velocity of light c. Thus, during accretion, at some point, mat- 
ter velocity should exactly match with the sound speed. This point is called a critical point 
or a sonic point. When matter crosses a sonic point, it becomes transonic. As Chakrabarti 
(C89ab, C90ab) pointed out, depending on the initial parameters, a flow may have multiple 
sonic points and therefore, depending on whether the shock conditions are satisfied or not, 
a flow may or may not have a standing shock. 

For the sake of completeness, we carry out the sonic point analysis by solving the above 
Eqs. (1-3) using a method similar to that used in C89b. 

3.1 Sonic point conditions 

In the present analysis, we use MISStrcss prescription (C96a) for computing and Wxcp is 
obtained from Eq. (3). For the accretion flow, the entropy equation (Eq. 4) can be simplifled 
as. 



Wx<p = -anil. 



(3) 



3 SONIC POINT ANALYSIS 



J — 1 p dx dx p. 

and then H{= / ph) takes the form, 

H = Ax{ga^ + 71?^) 



d IdP -iPdp Q- 




(4) 




(5) 
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where, 7 is the adiabatic index, A = —au^ and g = ^j^. Here, is the angular velocity 
of the accreting matter at the radial distance x, n is the polytropic index {n — :^), In and 
come from the definition of the vertically averaged density and pressure (Matsumoto 
et al. 1984). 

In the present analysis, we use Q- = 0, i.e., the cooling is ignored. This would be strictly 
valid if the accretion rate is low, so that the loss of energy by bremsstrahlung cooling is 
insignificant compared to the rest mass energy. 

After some simple algebra and eliminating da/dr etc. we get from the governing equations 
Eqs. l(a-c) and Eq. (2) the following first order linear differential equation: 



dt') 
dx 



N 



(6) 



where, the numerator is. 



N 



7X 

i92o2(5x - 3) 



2ix - 11 



2Q;n^A(a'5 + 7i?') + 



(7 + 1)^' 
(7-1) 



x(7 — l){x — 1) 
and the denominator D is 

^_ 2a^^ (7 + 1)^?^ 



■yx{x — 1) 
- Aau^{a^g + 71?^) 



+ 



x^ 



(2^ - 1) - 



7't?2 



(7) 



(8) 



(7-1) (7-1) 

Both and D are algebraic equations which makes this model easily tractable. 

At the sonic point, both the numerator and the denominator must vanish simultaneously. 
For D = 0, one can get the expression for the Mach Number M{xc) at the sonic point and 
is given by. 



M{xc 



-rub 



ml - Amarric 



2ma 



where. 



nia = -Aan7'(7 - 1)(2^ - 1) - 7(7 + 1), 
mfe = 27 - 2Aa-;,g-f{-f - l){g - 1) 
rUc = Aa^g'^i'y - 1). 

In the weak viscosity limit, an — ^ and the Mach number at the sonic point M{xc) ~ 
for an — > a result obtained in C89b. 
Setting A^ = 0, we get an algebraic equation for the sound speed at the sonic point which 
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is given by, 

■P'i^ci ^c) — 



anA{g + -fM^f ^ an^(5a; - 3){^ + 7M2} ^ M^{bx - 3) 



7X 



2XAM{g + 7M2) 



2 X 



7x(x — 1) 
2angA(g + -fM^) + 



x('-f — l)(x — 1) 
(7 + 



(7-1) 



(9) 



We solve the above quadratic equation to obtain tlie sound speed at the sonic point. Das 
et al. (2001) suggested that depending on a given set of initial parameters accretion flow 
may have a maximum of four sonic points where one of the sonic points always lies inside the 
black hole horizon for non-dissipative accretion flow. In our present study, we also expect a 
similar result. 



3.2 Nature of the sonic points 

A black hole accretion is always transonic. Thus the originally subsonic matter definitely has 
to pass through the sonic point to become supersonic before entering into the black hole. 
Depending on the initial parameters, a flow may have multiple sonic points. Nature of sonic 
point depends on the value of velocity gradients at the sonic point. It is easy to show that 
^ assumes two values at the sonic point. One of them is vahd for the accretion flow and the 
other is vahd for the wind. If both the derivatives are real and of opposite signs, the sonic 
point is saddle type. When the derivatives are real and of the same sign, the sonic point is 
nodal type. When the derivative is complex, the sonic point is spiral type (or 0-type, for 
non dissipative flow). See, C90b for details of classiflcations. In order to form a standing 
shock, the flow must have more than one saddle type sonic points. 

In Figure la, we plot the variation of speciflc angular momentum (Ac) as a function of the 
logarithmic sonic point location (xc) for a given viscosity parameter (an = 0.1). Here differ- 
ent curves are drawn for different specific energies at the sonic points. The energies, from the 
top curve to the bottom, are given by: Sc = 0.0007, 0.001, 0.003, 0.005, 0.007, 0.011, 0.015, 0.019, 0.023 
and 0.027 respectively. The long-dashed curve at the top represents the Keplerian angular 
momentum distribution which is completely independent of the initial flow parameters and 
depends only on the geometry. Sohd part of the curves represents the saddle type sonic 
points, dotted part of the curves represents the nodal type sonic points and the short- 
dashed part of the curves are for the spiral type sonic points. First notice that the sonic 
points always occur at angular momentum below Keplerian value. Notice that for lower 
values of specific energy at the sonic point, an accretion fiow contains all the three types of 
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log(xJ 

Fig. la: Variation of specific angular momentum (Ac) as a function of the logarithmic sonic point location (xc) for the viscosity 

parameter (an = 0.1). Long-dashed curve is the Keplerian angular momentum distribution. Solid curves represent the saddle 
type sonic points, dotted curves represent the nodal type sonic points and the short-dashed curves are for the spiral type sonic 
points. Shaded area is the nodal type sonic point region. 

sonic points in a systematic order: saddle — nodal — spiral — nodal — saddle for mono- 
tonic increase of location of sonic points. With the increase of energy Sc the region of spiral 
type sonic points gradually decreases and finally replaced by the nodal type sonic points 
though multiple sonic points still exist. Shaded area separates the nodal type sonic point 
region in the Ac — Xc plane. With further increase of energy all the nodal type sonic points 
also disappear and are replaced by saddle type sonic points. In this case, the flow has only 
one sonic point for a given sub-Keplerian angular momentum. Thus, for a given angular 
momentum of the flow, there exists a range of energy Smin < Sc < Smax such that the 
flow has multiple sonic points. In the inset, we zoom a small portion of the curve close the 
Keplerian value to highlight the fact that the angular momentum at the sonic point always 
remains sub-Keplerian when cooling process is ignored. In future (Das & Chakrabarti 2003, 
in preparation), we shall show that a flow can also be super-Keplerian when cooling is added. 

In Fig. lb, we show a very important aspect of a viscous transonic flow. Here we show how 
the angular momentum at the sonic point varies when viscosity parameter an is increased. 
We hold the energy at the sonic point to be fixed at Sin = 0.006. In the absence of viscosity 
(an = 0, the uppermost curve), the flow has all the three types of sonic points. Similar to 
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log(xJ 

Fig. lb: Variation of angular momentum at the sonic point as the viscosity parameter is varied. Specific energy at the inner 
sonic point is held fixed at 0.006. From the uppermost to the lowermost curve: an = 0,0-2,0.4,0.6 and 0.7 respectively. Other 

notations are the same as in Fig. la. 

Fig. la, here also we indicated the saddle, nodal and spiral type sonic points by the solid, 
dotted and short-dashed curves respectively. The uppermost long dashed curve represents 
the Keplerian angular momentum distribution. With the increase of an, more and more 
inner saddle type sonic points are replaced by nodal type sonic points and similarly nodal 
type are also replaced by spiral type sonic points. The curves, from the uppermost one to the 
lowermost one, are, for an = 0, 0.2, 0.4, 0.6, 0.7 respectively. For an = 0.7, all the inner 
saddle type sonic points disappear and only the spiral type points remain. Thus there exists 
a critical viscosity parameter an(c,j) at a given Sin for which all the inner saddle type sonic 
points are completely converted into spiral type ones. In this case, the flow has no choice 
but to pass through the outer sonic point only. Existence of such critical viscosities has been 
predicted in C90ab and C96a — below we compute their values exactly as a function of the 
inflow parameters. This behaviour also hints at the conclusion that the parameter space for 
the existance of a transonic flow may shrink with the increase of viscosity. 

We continue our investigation of the transonic nature of the flow and replotted a Figure 
similar to Fig. la but we increase an gradually while keeping the specific angular momentum 
at the inner sonic point to be fixed (Ac = 1.65). Values of an are, from the top to the bottom 
curve, an = 0, 0.25, 0.35, 0.5, 0.6 and 0.7 respectively. Sohd, dotted and short-dashed lines 
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Fig. Ic: Variation of the specific energy at the inner sonic point as a function of the viscosity parameter an- Prom the 
uppermost curve to the lowermost curve, an = 0, 0.25, 0.35, 0.5, 0.6 and 0.7 respectively. Other notations are the same as in 

Fig. la. 

represent the saddle type, nodal type and the spiral type sonic points respectively. Long- 
dashed line separates the positive and negative energy regions in the Sc — Xc plane. Notice 
that, for increasing ctn, saddle type sonic points are gradually replaced by the nodal and 
spiral type sonic points: outer saddle type sonic points recede further away and the inner 
saddle sonic points proceed toward the black hole horizon. For an = 0.7, the inner saddle 
type sonic points completely disappear and become spiral type. This behaviour points to 
a critical value of viscosity parameter (ctn.c) which separates the accretion flow from the 
multiple sonic point regime to the single sonic point regime at a given Ac. It is also clear 
that at the same sonic point, specific energy steadily decreases for increasing an- This is 
because, when an is increased, the accreting matter tends to become a Keplerian disk closer 
to the black hole and becomes more strongly bound with lower energy. Note that the energy 
at the outer sonic point remains always positive for the all initial parameters. 

In our final study of the nature of the sonic points, we chose n, the polytropic index, to 
be our free parameter. For a very relativistic flow, or a radiation dominated flow, n — 3, 
but for a mono-atomic, non-relativistic gas, n — 3/2. In Fig. Id, we show the variation of Ac 
with sonic point location Xc- We keep the specific energy at the sonic point to be 0.001. The 
long-dashed curve is the Keplerian distribution as before. We note that with the increase of 
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Fig. Id: Variation of the specific angular momentum at the inner sonic point as a function of the polytropic index n (marked 
on each curve). Generally, the number of sonic points decreases by decreasing n. Specific energy has been kept fixed at 0.001 
except for the dashed curve where it is 0.0015 to show that for a given polytropic index, number of sonic points increases with 
decreasing energy. 

the adiabatic index 7, i.e., decrease of the polytropic index n, the number of sonic points 
decreases from three to one. In this example, in the extreme non-relativistic regime (n = 1.5) 
the accretion flow has a single saddle type sonic point for any specific angular momentum. 
For the same energy, for n — 1.75 there are three sonic points, indicating that a standing or 
an oscillating shock in the flow may be possible. In this Figure, we also show that for a given 
n (such as for n — 1.75) if we increase energy at the sonic point {£c — 0.0015), multiple sonic 
points disappear and a single sonic point forms. This indicates that there must be a critical 
value of Ec associated with each n above which multiple sonic points do not exist when all 
other parameters arc kept fixed. 

The general behaviour of the fiow at the sonic point is best seen in Fig. 2, where we 
depict the surface F{Sc, Ac, Xc) — (Eq. 9) for ctTr = 0.01. Sonic points Xc are plotted along 
X-axis in the logarithmic scale, Ac is plotted along the Y-axis and Ec is plotted along the 
Z-axis. At high angular momenta, there are three sonic points, but they merge to become 
one at lower angular momenta. Below a critical value Ac, the flow does not have more than 
one sonic point. 
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Fig. 2: Gradual change in the number of physical sonic points is easily seen in this three-dimensional view of F{£c, Xc,Xc) = 
(Eq. 9) surface. At high angular momenta, there are three sonic points, but they merge to become one at lower angular 
momenta, an = 0-01 has been chosen. 

4 GLOBAL SOLUTION TOPOLOGY 

Basic criteria for studying shock properties is that the accretion flow must have a multiple 
saddle type sonic points and the shock should join two solutions — one passing through 
the outer sonic point and the other passing through the inner sonic point. The solution 
topologies have been already discussed in C96a. The current paper studies in great detail 
the topologies associated with / = 1 of C96a. In particular, we show below new pathways 
through which topologies may vary. 

In Fig. 3(a-b), we have shown how the flow topologies change with the viscosity param- 
eter ctn and the specific angular momentum Aj^ at the inner sonic point Xin- Two distinct 
behaviours have been highlighted here: one at a low angular momentum regime (Fig. 3a) 
and other at a high angular momentum regime (Fig. 3b). In Fig. 3a, we keep the inner sonic 
point fixed at Xin = 2.665 and the specific angular momentum at this point is Aj„ = 1.68. 
At low angular momentum and without viscosity (the box at extreme left in Fig. 3a) the 
subsonic flow enters into the black hole after passing through the inner sonic point. In the 
second box, viscosity is slightly higher and topologies are closed for the same inner sonic 
point. So, for a given set of parameter, there must be a critical viscosity parameter (anc) for 
which open topologies become closed ones. We will discuss critical viscosity rigorously in the 
§7. Accretion with parameters causing this kind of topology never joins with any Keplerian 
disk unless a shock is formed (This will be shown below.). When a standing shock formation 
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Fig. 3(a-b): Variation of global solution topologies of the viscous a<;cretion flow around black holes. In (a), drawn for \{xin = 
2.665) = 1.68, four panels show how open topology at lower viscosity becomes closed at higher viscosity. In (b), drawn for 
-^(^in = 2.359) = 1.78, closed topology opens up again. 

is not possible, an accretion flow passes through the outer sonic point directly before falling 
in to the black hole. For a further increase of an (next two boxes) closed topology shrinks 
gradually and finally disappears leaving behind only the outer sonic point (Bondi Type). 
This is directly analogous to the shrinking of the phase space of a simple harmonic oscillator 
in presence of damping (C90a). These solutions are basically the same as the f — 1 case of 
Fig. 2a of C96a. 

In Fig. 3b, where solutions arc plotted with a higher specific angular momentum at the 
inner sonic point {xin — 2.359 and Ai„ = 1-78), explanations of first and second box is similar 
to earlier one (Fig. 3a), but in the third box (ctn = 0.05), accretion flow topology reverts 
its direction of spiraling and the flow can join with a Keplerian disk very close to a black 
hole. All the differences between this two flgures (Fig. 3a. and Fig. 3b) is mainly due to 
the difference of specific angular momentum at the sonic point rather than the change of 
sonic point locations. In Fig. 11 (below), we will show that the nature of the accretion ffow 
topologies have strong dependence on angular momentum at the sonic point. For a higher ctn 
(next two boxes) Keplerian disk comes even closer to the black hole and topologies passing 
through the inner sonic point becomes Bondi type. We suspect that two limits of viscosity 
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Fig. 4: Variation of global solution topologies of the viscous accretion flow when the outer sonic point is kept fixed at Xout = 30rg 
while the viscosity parameter is varied. The closed topologies (with one saddle type and a spiral type at the center) at lower 
an become open. Eventually, the saddle type also disappears to produce a nodal type sonic point. 

parameters would cause an oscillation of the inner part of the Keplerian disk, but we cannot 
be certain about it without a time dependent numerical simulation. 

In order to show that flow topologies often can take new pathways than what was already 
known (C96a), we draw in Fig. 4, we have plotted solution topologies passing through the 
outer sonic point chosen at Xout = '^5rg and specific angular momentum \{xout) is 1.8. The 
viscosity parameter is varied (marked on each box). For a lower a^, the topologies are closed 
as in Fig. 4 of C96a and the flow having this topology cannot be transonic anyway. When an 
is increased, closed topologies gradually open up (unlike C96a where Xout — 35 was chosen 
and the opening of topologies did not occur) and if shock condition is satisfied, the accretion 
flow passing through the outer sonic point jumps in to the subsonic branch and go through 
the inner sonic point before entering into the black hole. For higher an, the same outer sonic 
point no longer remains saddle type. First it becomes of nodal type and then it becomes of 
spiral type (Fig la). Considering that the outer sonic point recedes farther away with the 
increase of an, this behaviour is not surprising. This, together with Fig. 4 of C96a show that 
there could be more than one way of reaching nodal topology. 
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Fig. 5: Classification of the parameter space spanned by the specific energy and the specific angular momentum of the flow at the 
inner sonic point. The bounded regions drawn for different viscosity parameters (marked at the inset) contain allowed solutions 
which may pass through the inner sonic point. As viscosity increases, the region shifts towards lower angular momentum and 
higher energy. 

5 CLASSIFICATION OF THE PARAMETER SPACE 

An important part of understanding a viscous flow is to classify the parameter space as a 
function of the viscosity parameter. In Fig. 5, we have separated the parameter space for 
the accretion flow which can pass through the inner sonic point. Angular momentum at this 
sonic point (Aj„) is varied along the X-axis and the corresponding speciflc energy at the inner 
sonic point (Sin) is plotted along the Y-axis. The region bounded by a given curve contains 
parameter space in which multiple sonic points are possible. For instance, for an = 0, the 
region bounded by the solid curve is identical to the region found by C89b. For increasing 
ctn, the region of multiple saddle type sonic point is reduced near high angular momentum 
side while it increases in the lower angular momentum side. It may be recalled (C89b and 
Fig. la) that at low angular momentum, the number of sonic point is just one. With the rise 
of an, the angular momentum at the sonic point is increased, increasing the number of sonic 
point. At a higher angular momentum, the situation is just the opposite. In this case, there 
are already multiple sonic points for an ~ and for high enough an, viscosity transports 
angular momentum very rapidly causing a steep rise in angular momentum itself. This, in 
turn, means that the flow can have only one saddle type sonic point in this case. 
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6 STANDING SHOCKS AND FURTHER CLASSIFICATION OF THE 
PARAMETER SPACE 

Though in general astrophysical contexts, shocks are ubiquitous and possibly non-stationary, 
in an accretion flow, the location, strength and thermodynamic quantities could be quantifled 
very exactly by using the Rankine-Hugoniot conditions (RHCs). The study is similar to the 
study of shocks in solar winds (e.g., Holzer & Axford, 1970) and white-dwarf surfaces. 

The shock conditions which we employ here are the usual RHCs presented in C89b, i.e., 
(a) the local energy flux is continuous across the shock; (b) the mass flux is continuous 
across the shock (c) the momentum balance condition is satisfled and flnally, (d) angular 
momentum should be continuous across the axisymmetric shock. 

The way an accretion flow moves around a black hole, as seen from local rotating frame, 
is as follows. First the flow, subsonic at a very far away, passes through the outer sonic point 
and becomes supersonic. The RHCs then decide whether a shock will be formed or not. 
Of course, our consideration of satisfying RHCs at a given location holds only if the shock 
is thin, i.e., viscosity is low. Nevertheless, we continue to use this prescription at higher 
viscosities to have a flrst order guess of the shock location. Similarly, we assume that there 
is no excess source of torque at the shock itself, so that the angular momentum may be 
assumed to be continuous across it. This condition may be violated when magnetic flelds 
are present. In the presence of large scale poloidal magnetic flelds, there could be magnetic 
torques which could make flow angular momentum discontinuous. 

6.1 Method of calculating the shock locations 

Accretion flow after passing through outer sonic point jumps (shock) into subsonic branch 
and becomes supersonic while crossing inner sonic point before falling into the black hole. In 
our present study we begin numerical integration from inner sonic point and proceed towards 
the outer edge of the accretion disk to look for shock location. During integration along the 
subsonic branch, it is possible to calculate all the local variables {i.e.-d, a, M, p at the post- 
shock region, in terms of the initial flow parameters. We calculate total pressure, local flow 
energy, speciflc angular momentum at shock using these subsonic local variables. At the 
shock, total pressure, local flow energy, mass accretion rate (one of the flow parameters) and 
speciflc angular momentum are conserved. These conserved quantities at the shock give the 
other set of supersonic local variables for the supersonic branch. This supersonic set of local 
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variables help to get outer sonic point uniquely for a accretion flow with fixed inner sonic 
point and other initial flow parameters when integration takes place towards the outer edge 
of the black hole. Thus, the accretion flow can be connected with both the saddle type sonic 
points through shock for dissipative system and this determines the standing shock location 
for a given set of initial parameters. 

We compute the supersonic local flow variables in terms of subsonic local flow variables 
in the following way: 

Our model accretion flow is in vertical equilibrium and total pressure of the accretion 
flow at any given point is given by, 

^ = ^y + E^?^ (10) 

where, W and S are the vertically averaged thermal pressure and density respectively. 

We use the mass conservation equation (Eq. lb) in the Eq. (10) and calculate the sound 
speed a in terms of radial velocity (■»?) at the shock Xs in the supersonic branch which is 
given by, 

a'' =CiC2^-C2^^, (11) 

where, Ci = and 62^^. 

Prom the local flow energy equation, radial velocity at shock in the supersonic branch 
can be calculated using Eq. (11) and is given by. 



where, = 2nCiC2, A = -2nCiC2, d^ = 1E- + and E is the local flow energy. 

Here, total pressure and flow energy at the shock is calculated with the help of the 
subsonic flow variables. We consider only the '+' sign as we interested to get the local flow 
variables in the supersonic branch. This radial velocity is used to get the sound speed (Eq. 
11) in the supersonic branch of the flow. Now we use supersonic flow variables (■»? and a) to 
get outer sonic point by numerical integration and it completes the accretion flow solution 
having shock. 

In Fig. 6a, we have shown how shock may be formed by joining two solutions, one with 
a lower entropy passing through the outer sonic point, and the other with a higher entropy 
passing through the inner sonic point. Two curves are drawn for two different energies 
(marked). Sin and E^t are the energies at the inner and outer sonic points for a shocked 
accretion flow which has a standing shock. Due to viscous heating processes, energy is 
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Fig. 6a: An example of how a standing shock might form in a viscous transonic flow is depicted here. Flow passing through the 
outer sonic point at Xout = 39.7 and energy £out = 0.00608 has a shock and passes through the inner sonic point at Xin = 2.78 

where its energy is £i„ = 0.007555. Inset shows the details. 

increased and the shock wave is formed when flow jump from the lower energy solution to 
higher energy solution. If we included only the cooling process, the situation would have 
been reversed. The flow parameters are Xout — 39.7, Aj„ = 1.65, an — 0.05 and 7 = 4/3. The 
shock condition uniquely determines the inner sonic point, which is at Xj„ = 2.78. The end 
positions of the long arrow mark the locations of the sonic points. In the inset, we zoom a 
selected region in the Ac — Xc plane to show explicitly that the angular momentum is indeed 
decreased. 

In Fig. 6b, we present the complete solution of the flow which includes a standing shock 
in a viscous flow for the same set of parameters used to draw Fig. 6a. Arrows indicate the 
direction of the accreting flow. Subsonic accreting flow passes through the outer sonic point 
(O) and becomes supersonic. At Xg, shock conditions are satisfied — the flow jumps from 
the supersonic branch to the subsonic branch and subsequently pass through the inner sonic 
point (1). In this particular case, the shock conditions are satisfied at Xs = 18.2 and shock 
is denoted by the dotted vertical line. 
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Fig. 6b: Actual solution topology for the case discussed in Fig. 6a is shown here which, along with the outer (O) and inner 
(I) sonic points, also showed the shock transition at Xs = 18.2 (vertical dotted line). The arrowed curve is followed by a flow 
while entering into the black hole. 




Fig. 7: Variation of the region of the parameter space which forms a standing shock as a function of the viscosity parameter 
an . The region shrinks with the increase of viscosity parameter. 
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Fig. 8: Mapping of tlie parameter space of the pre-shock region (solid boundary) spanned by (fin, Afn) onto the parameter 
space in the post-shock region (dashed boundary) spanned by {£out,^out) in a viscous flow (an = 0.01). 

6.2 Parameter space which allows standing shocks 

In Fig. 5, we have already classified the parameter space in terms of the number of sonic 
points in the flow. Presently, in Fig. 7, we concentrate on the region which allows only 
standing shocks in a viscous flow. The viscosity parameters are marked. The region marked 
au — coincides with that in C89b and in Chakrabarti (1996b, 1998; hereafter C96b and 
C98 respectively) when appropriate models are considered. Compared to the inviscid case, 
the effective region of the parameter space shrinks in the high angular momentum side 
when the viscosity is increased. The situation is exactly the opposite at the lower angular 
momentum side. We observe that even at angular momentum as low as 1.4, standing shocks 
could be formed if the viscosity is high enough. Above a critical viscosity (which depends 
on other parameters as will be shown in Fig. 13 below), this region disappears completely. 

We continue our study of the parameter space which may allow multiple sonic points. In 
Fig. 8, we show a curious feature: the mapping between the post-shock parameters and the 
pre-shock parameters. We plot the region of the post-shock parameters at the inner sonic 
point {Sin, Kn) (bounded by the solid curve) and the region of the pre-shock parameters 
at the outer sonic point {Sout, ^out) (bounded by the long dashed-curve) for a shock which 
is determined through RHCs for au — 0.01. For each and every point in the pre-shock 
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parameter space region, there exists a point in the post-shock parameter space region and 
therefore we have a complete solution. For definitiveness, we also show vertical dashed and 
dotted lines in the two different angular momentum range in the post-shock parameter space 
region in which Aj„ is kept fixed but Sin is varied. The corresponding pre-shock parameters 
form a curve, indicating that both the angular momentum and the energy had to be adjusted 
to get the self-consistent solution. In a non-dissipative flow, there is no variation of energy 
and angular momentum in the accretion flow. As a result, both the inner and the outer sonic 
point parameter spaces merge (C89b, C98). 

6.3 Parameter space which may allow oscillating shocks 

An important role that is played by oscillating shocks is to produce the so-called quasi- 
periodic variations of X-ray intensity from galactic black hole candidates. In the inviscid flow, 
the region of the parameter space which produced multiple sonic points, and yet, RHCs were 
not satisfied was important for such type of oscillating shocks (Ryu, Chakrabarti & Molteni, 
1997). Here, the winds are also produced sporadically from the post-shock region. In presence 
of coohng, especially when the cooling time scale roughly agree with the infall time-scale, 
otherwise standing shocks the shocks may oscillate (Molteni, Sponholz and Chakrabarti, 
1996). 

It is to be noted that it is, in a general flow, very difficult to divide the parameter space 
in terms of whether the shock will exist or not. This is because, when there are shocks, 
at least the RHCs allow us to map the pre-shock and the post-shock flow parameters (see. 
Fig. 8). But when there are no shocks, it is not straight forward to map these two sets of 
parameters. Thus, one has to rely on global topological behaviour of the flow solutions and 
whether they allow multiple sonic points. 

One of the criteria to use is to check which parameter space allows one to have the inner- 
sonic point energy to be larger compared with the outer sonic point energy. For instance, 
in Fig. 9a, we plotted the inner sonic point energy (Sin) along X-axis and outer sonic point 
energy (Sout) along the Y-axis for a set of inner sonic point angular momentum (Xin) when 
accretion flows pass through shocks. It is clear that Sout varies almost linearly with Sin and 
the nature of this variation depends only on Xin. It is not unwarranted to assume that a 
similar linear variation will continue for shock-free solutions also at least if the viscosity is 
low. Thus we extrapolate this variation for the shock-free solution for the lower values of 
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Fig. 9a: Example of variation of the outer sonic point energy Sout as a function of the inner sonic point energy Sin when the 
flow has a shock, an = 0-1 and X{xin) = 1.50 for the topmost curve. Curves have an increment of AX{xi„) = 0.02 while going 

towards the bottom. 

Sin till when Sout ^ for the same Ai„ keeping in mind that accretion flow topology passing 
through the inner sonic point must remain closed (Fig. 3a-b). We follow this procedure to 
estimate the cut-off Sin for different Aj„ and obtain the region in the parameter space where 
accretion flow has more than one X-type sonic point. 

In Fig. 9b, we show the division of the parameter space {Sin, Kn) for the viscosity pa- 
rameter an = 0.01 on the basis of number of sonic points. Flows with parameters from the 
region ACD have more than one X-type (saddle type) sonic points and RHCs are also satis- 
fied. Flows with parameters from the region ABC have more than one X-type sonic points 
but RHCs are not satisfied here. From our previous experience with non-dissipating fiows, 
we predict that those solutions with multiple sonic point which do not produce standing 
shocks must be producing oscillating shocks. This region becomes bigger when the viscos- 
ity parameter is reduced. Rest of the parameter space gives solutions with closed topology 
passing through the inner sonic point. 

6.4 Parameter space for all possible solutions 

Figure 10 shows the classification of the parameter space in the energy-angular momentum 
{Sin, Kn) plane in terms different accretion fiow topologies (small box) for an = 0.01. Solid 
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Fig. 9b: Division of the parameter space (Sim^in) for a viscosity parameter (an = 0.01) on the basis of number of sonic 
points. Region separated by the dotted line has more than one X-typc (saddle type) sonic points and flows in this region form 
standing shocks. The region surrounded by the dashed courve has more than one X-type sonic points but the Rankine-Hugoniot 
conditions are not satisfied here. 




Fig. 10: Division of the parameter space (£in,^in) (marked) for a viscosity parameter (an = 0.01) on the basis of solution 
topologies shown in boxes (marked). See text for details. 
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boundary separates the region in the parameter space for closed topologies passing through 
the inner sonic point in general. Further sub-classes are indicated by dotted, dashed and dot- 
dashed curve which classified the solution topologies depending on their behaviour. Examples 
of solution topologies with initial parameters taken from different regions (marked) of the 
parameter space are plotted in seven small boxes (marked) . All the small boxes depict Mach 
number variation as a function of the logarithmic radial distance. The box labeled S shows 
an accretion flow solution which passes through a shock. Dotted vertical line with an arrow 
indicates the location of the standing shock. The solution drawn in the box marked OS 
is an accretion fiow which has multiple sonic points but does not satisfy RHCs after the 
flow becomes supersonic. Prom our earlier experience with an inviscid flow, this topology 
is expected to give rise to an oscillating shock solution. The box marked OAC shows a 
new type of solution topology having multiple sonic points. One branch of the topology is 
closed and the other branch is open. This kind of solution is available in a small region of 
the parameter space shown in the inset on the upper-right corner. Solutions inside the CIl 
box have closed topology (inner spiral going anti-clockwise) having only one saddle type 
sonic point and this kind of solution belongs to a large region of the parameter space with 
relatively lower angular momentum region. The box CI2 shows a similar result as CIl but 
here the nature of the topology is different (inner spiral going clockwise) and this type of 
solution exists at higher angular momentum region. The box labeled I* represents accretion 
flow solution which only passes through the inner sonic point. This solution could be for an 
accretion or wind and the initial parameters for this type of topology belongs to the region 
indicated by I* in the parameter space. Topology with parameters taken from O* region 
of the parameter space is also plotted in the box marked O*. An accretion flow solution 
with these parameters passes only through the outer sonic before falling into the black hole 
(similar to a Bondi flow). 

Since the strength of the shock determines the jump in temperature and density, it may 
be worthwhile to study the shock strength. We deflne this as the ratio of the pre-shock Mach 
number to the post-shock Mach number. As an example, in Fig. 11, we show the variation 
of the shock strength as a function of viscosity parameter an. This Figure is drawn for 
Xin — 2.795 and Ai„ = 1.65. For lower viscosity limit, the strength of the shock is weak. It 
increases smoothly with the gradual increase of viscosity and there is a cut off at a critical 
viscosity limit where the shock disappears. 
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Fig. 11: Variation of the ratio of the pre-shock to post-shock Mach numbers as a function of viscosity parameters for a fixed 
set of initial conditions {xi„ = 2.795 and Ai„ = 1.65). The shock disappears beyond the critical parameter an ~ 0.017. 

7 DEPENDENCE OF THE CRITICAL VISCOSITY PARAMETER 

In our earlier discussion, we already hinted that there must be a critical viscosity parameter 
for which the flow topology must change its nature from an open topology to a closed 
topology. Presently, we quantify the nature of this critical viscosity. Indeed, we find that 
there are in effect two critical viscosity parameters: one at the boundary which separates 
the closed topology from the open topology while the other splits region of closed topology 
in terms of whether shocks can form or not. Not surprisingly, these are infiow parameter 
dependent, and thus do not have universal values. Nevertheless, these are useful, since they 
give us insights into the cases in which shocks may be possible. 

Figure 12 shows the variation of the critical viscosity parameters with the angular mo- 
mentum at the inner sonic point. Different regions are marked. We note that for higher 
viscosity parameters (a^), shocks are formed in the lower angular momentum domain. As 
the angular momentum is increased, the shock disappears. This was also expected from our 
discussion of Fig. Ic. 
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Fig. 12: Critical viscosities separating standing from oscillating shocks and closed topologies from open topologies. 
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Fig. 13 (a-b) : Variation of shock location with viscosity parameter and (a) inner sonic point and (b) specific angular momentum. 
Shock location always decreases with increase of viscosity till the critical viscosity parameter is reached beyond which the shock 
ceases to exist. 
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8 DEPENDENCE OF THE SHOCK LOCATION ON VISCOSITY 
PARAMETER 

In our study of shock properties, we have aheady mentioned that the shocks disappear when 
viscosity is more than a critical value. In Fig. 13 (a-b), we show how the shock location 
depends on the viscosity parameter when the other two free parameters, i.e., the shock 
location {xin) and angular momentum (Aj„) are kept fixed. In Fig. 13a, variation with the 
inner sonic point is shown when the angular momentum is kept fixed, while in Fig. 13b, the 
variation with the angular momentum is shown keeping the inner sonic point fixed. In all the 
cases, the shock location is reduced with the increase in viscosity parameter till the critical 
viscosity parameter is reached beyond which the shock ceases to exist. This is significant 
because in an accretion flow, when viscosity is increased, the accretion rate is also increased 
and a black hole candidate goes from spectrally hard to spectrally soft state (CT95). Thus 
if the shock oscillation is indeed the cause of quasi-periodic oscillations (QPOs), then the 
frequency should increase with the accretion rate and finally as the shock ceases to exist, the 
QPOs also should disappear in softer states. Observation of such features could be used to 
verify if the shock oscillations may be the prime cause of the QPOs in black hole candidates. 

9 CONCLUDING REMARKS 

In this paper, we have extended our earlier results of the study of shock formation to include 
a very difficult yet more realistic case of viscous polytropic flows. Some of the results have 
been touched upon in C96a but the new results in our work include a detailed study of 
the parameter space in which shocks form even in presence of viscosity. We found a large 
number of important results: (a) That there exists two critical viscosity parameters which 
separate the region of the parameter space in three parts - 1) in which the flow has a Bondi- 
type single sonic point; 2) in which there are three sonic points but no Rankine-Hugoniot 
relations are satisfied and 3) when Rankine-Hugoniot relations are satisfied. These critical 
viscosity parameters decrease with the increase of the specific angular momentum of the 
flow at the inner sonic point, (b) That at high viscosities, standing and oscillating shocks 
may form if the flow has very little angular momentum at the inner sonic point, while at 
low viscosities the situation is exactly the opposite. It is widely believed that accreting 
matter on galactic and extra-galactic black holes could be of very low angular momentum, 
especially when the central compact object is accreting winds from the nearby star or stars. 
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This brings out the possibiUty that shocks may be active ingredients of an accretion flow. 
Our results, with a very plausible accretion flow models, indicate that the standing and 
oscillating shocks are produced even for large viscosity parameters, (c) That the shock 
location is reduced with enhancement of viscosity parameter. This, coupled to earlier results 
(Chakrabarti & Manickam, 2000) that the infall time is proportional to the period of quasi- 
periodic oscillation (QPO) of X-rays from black holes imply that the QPO frequency should 
increase as the viscosity is increased. This is consistent with the observational findings that 
QPO frequency is increased as the spectral slope softens, widely known to be due to increase 
in viscosity and accretion rate. 

One of the questions we have not addressed here is the stabihty properties of these shocks. 
A number of authors pointed out that while the shocks are stable, they should undergo oscil- 
lations, either radially, or vertically, or non-axisymmetrically (Molteni, Toth & Kuznetsov, 
1999; Gu & Foglizzo, 2003). We anticipate that our shock solutions in viscous fiows would 
suffer similar types of oscillations, especially when the viscosity is low. In particular, Gu & 
Foglizzo (2003) while studying shocks in inviscid, isothermal flows found such instability and 
interpreted as due to cycles of acoustic waves between their co-rotation radius and the shock. 
In their interpretation this could be a form of Papaloizou-Pringle instability (Papaloizou & 
Pringle, 1984) which is known to destabilize accretion tori when the angular momentum gra- 
dient is less than a certain value. If so, such an instability could disappear at high enough 
viscosity. This could have a bearing on the quasi-periodic oscillations of observed X-rays in 
galactic and extra-galactic black hole candidates in that QPOs would cease to exist above 
a certain frequency. The interesting aspect is that these so-called 'instability' only causes 
oscillation of shocks and does not destroy the shock (Molteni, Toth & Kuznetsov, 1999). 

This work is partly supported by a project (Grant No. SP/S2/K-15/2001) funded by 
Department of Science and Technology (DST), Govt, of India. 
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